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ABSTRACT 

There have been many proposed explanations for the larger-than-expected radii of some tran- 
siting hot Jupiters, including either stellar or orbital energy deposition deep in the atmosphere 
or deep in the interior. In this paper, we explore the important influences on hot- Jupiter ra- 
dius evolution of (i) additional heat sources in the high atmosphere, the deep atmosphere, and 
deep in the convective interior; (ii) consistent cooling of the deep interior through the planetary 
dayside, nightside, and poles; (Hi) the degree of heat redistribution to the nightside; and (iv) 
the presence of an upper atmosphere absorber inferred to produce anomalously hot upper at- 
mospheres and inversions in some close-in giant planets. In particular, we compare the radius 
expansion effects of atmospheric and deep-interior heating at the same power levels and derive 
the power required to achieve a given radius increase when night-side cooling is incorporated. 
We find that models that include consistent day/night cooling are more similar to isotropically 
irradiated models when there is more heat redistributed from the dayside to the nightside. In 
addition, we consider the efficacy of ohmic heating in the atmosphere and/or convective interior 
in inflating hot Jupiters. Among our conclusions are that (i) the most highly irradiated planets 
cannot stably have uB > 10 km s _1 • Gauss over a large fraction of their daysides, where u is the 
zonal wind speed and B is the dipolar magnetic field strength in the atmosphere, and (ii) that 
ohmic heating cannot in and of itself lead to a runaway in planet radius. 

Subject headings: planetary systems - radiative transfer 



1. INTRODUCTION 

The discovery of the "hot Jupiter" class of exo - 
planets (jMavor fe Quelozlll995l:[Marcv fe Butlerlll996[) 
ushered in a variety of mysteries, of which perhaps 
the first was how they came to exist, since planet- 
formation models did not gener ally predict Jupiter- 
mass objects in few-day orbits^ (|Bossl H995t iLissauerj 
Il995f) . Whe n HP 209458 b the fi rst transiting planet 
disco vered (jHenrv et al.l I2000t ICharbonneau et"aLl 
2000), was found to have a radius <; 30% larger than 
Jupiter's, a new mystery was born: why are some 
hot Jupiters as inflated as they are (iBurrows et ail 
2000 ; iFortnev fc Hubbardl l200i ILauehlin et al.l 120051: 
Burrows et alj l2007aD ? More than a decade later, 



1.0 Rj, where Rj is Jupiter's radius0 The largest of 
the hot Jupiters have stri kingly large radii, incl uding 
the 2.04-i?j HAT-P -32b (IHartman et all I2011D. the 



there were at least 50 transiting planets known with 
masses greater than a third of Jupiter's and radii 
greater than 1.3 times Jupiter'sB The radii of these 
objects were initially seen as surprising because the- 
oretical evolutionary models suggested that, at the 
inferred ages of the systems (often more than a bil- 
lion years old), a hydrogen/helium gas-giant planet 
"should" have a radius significantly smaller, close to 



dave@ias.edu, burrows@astro.princeton.edu, 
1 Although, Struvc ( 1952) did suggest the possibility that 
objects such as what are now called "hot Jupiters" might exist. 
See http://cxoplanct.cu (Schneider et al. 2011), 



http://exoplancts.org (Wright et al. 2011), 
or http: / /cxoplanct.hanno-rein. dc/ (Rein 2012). 



1.99-Rj WASP-1 7b (lAnderson et al.ll201lD . the 1.83- 
Rj HAT-P- 33b (IHartman et all 1201 lh. the 1.74-Rj 
W ASP-12b (IChan et al.ll2011h . and the 1.71-Rj TrES- 
4 (jChan et al.l 1201 lh . At present, although a number 
of potential explanations for these planets' large radii 
have been suggested, it is unclear whether a single "in- 
flating" mechanism predominates. 

The processes that have been suggested to ex- 
plain the surprisingly large radii include both (i) ex- 
tra power sources in the interior and (ii) enhanced 
opacity or atmospheric s tratification (|Burrows et al.l 
2007a; Baraffc et al. 2010). In addition, there is a sub- 
tle distinction between the radius typically calculated 
in evolutionary cooling models (the radius of the pho- 
tosphere) and the radius observed in tr ansit; however, 
this s o-called "transit ra dius effect" (jBurrows et al.l 
2003; Baraffc et al. 2003) increases the apparent ra- 
dius by no more than ~5% relative to the radial pho- 
tosphere. Explanations of both category (i) and cat- 
egory (ii) tend to rely on the influence of the star 
and, therefore, naturally predict that the hottest hot 
Jupiters (with incident irradiation more than 10,000 

3 We note that some planets are smaller than expected, which 
can be explained by the presence of either a heavy-element core 
or heavy elements thr oughout the envelope dGu illot et al. 2006; 
IBurrows et~aT]|2007all . 
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times Jupiter's) should have the highest probability 
of being significantly larger than Jupiter. This predic- 
tion is borne out by obse rvations ([Laughlin et al.l2011t 
iDemorv fe Seagerll201lD . as can be seen in Fig. Q] 

Heating the convective interior of a planet di- 
rectly in creases its bu lk entr op y and, hence, its 

Bodcnhe imer et al.1 



radius (Burrows et al 



Il99i 



200H: ILiu et al.l 120081: ISpiegel fc Burrows! I201S.. 
Marleau fc Cumming||2013l) . Possible sources of heat 



in the interi or include tidal d i ssipat i on of orbit al 
eccentricity (iBodenheimer et all 120011: IWul l2005allbl 
Jackson et all 120081: ILiu et al.l l2008t llbeui fc Burrow: 



2009; 



iMiller et al.1 120091 : ILeconte et all I20T 
Ibgui et al.l 120101. 120111) . dissipation of thermal 
tides (|Arras fc SocratesT l2009bllal [20Tot) . dissipa- 
tion of downwardly propogating gravity wave s 
(jGuillot fc Showmaril2T)02l :ISh owman fc Ouillotif2 ()02). 
and ohmic heating from the dis sipation of currents in 
the partially ionized interior dBatvgin fc Stevenson! 
l20ia IBatvgin et afl [2MljBuang fc Cummind I2012t 



see also iPerna et al.l 2010a, who discussed a similar 
ohmic heating mechanism that might occur in the 
radiative portion of the atmosphere). Since many 
hot Jupiters have nearly circular orbits in which 
there should be little present-day tidal dissipation, 
steady-state tidal heating probably cannot generically 
explain all inflated radii. The possibility remains 
that a planet might have experienced an episode 
of tidal circularization in the recent past, injecting 
significant energy, increasing its radius, and effectively 
"resetting the clock" o n its evolutionary shrinkage 
(|Ibgui fc Burrows] 12009! ) . However, since the Kelvin- 
Hclmholtz timescale of an inflated hot Jupiter planet 
is on the order of ~100 million years, far shorter 
than the ages of the systems, this sort of explanation 
requires some fine tuning in the timing and, therefore, 
in the initial orbital configuration and efficiency of 
tidal dissipation. While it might explain the radii 
of some inflated hot Jupiters, it is n ot a preferred 
explanation for the ent ire popula tion (ILeconte et all 
120101 : llbgui eTaIll20ll . Indeed. Hbgui et all (|201lT 
found via an exhaustive parameter study that the 
extremely large radius (~1.7 Rj) of TrES-4, for 
instance, cannot be explained simply by a recent 
episode of tidal dissipation. 

Heating the atmosphere of a planet can puff up 
its outer few scale heights, but cannot directly ex- 
plain the radii of many of the larger hot Jupiters 
(e.g., those with radii >1.5 Rj must almost cer- 
tainly have inflated interiors). However, since the at- 
mosphere of a planet is the conduit through which 
the convective interior radiatively loses its heat (en- 
tropy), the thermal structure of the atmosphere crit- 
ically mediates the evolution of the thermodynamics 
of the interior. The intense stellar irradiation that 
hot Jupiters experience can dramatically change the 
vertical profiles of their atmospheres, inducing a deep 
isothermal layer through which evolutionary cool- 
ing is largely effectively stanched dGuillot et alJP L 996; 
Burrows et al.ll200a[20Ml2006l[2T^7al:lChabrier et al l 
2004: ILiu et alll2008t " THans en 2008; Guillotll2010h . and 



significantly modifying the evolution from that of a 
less strongly irradiated planet (s uch as those of our 
solar system; iFortnev et al.l 1201 ID or of widely sepa- 
rated planets or b r own dwarfs (Burrows et all [20011: 
i Baraffe et al.l 120031: iMarlev et all 120071: IFortnev et al l 
l2008bt iBurrows et al l 120111 : ISpiegel fc Burrows! 120121 



aJjrAJll; 

iMarleau fc CummirigTT2013l : fPaxton et al.l 12013! ) . The 
loss of heat from the interior can be further slowed 
by in creased atmospheric opacity ()Burrows et al.l 
2007a), or by the reduced efficiency of convective 
heat transport due to double-diffusive convectio n 
(jGhabrier fc Baraffdl27)07T : ILeconte fc Chabrierl[20Tl . 
Moreover, extra heating in the atmosphere of a planet, 
ultimately driven by the irradiation, can reduce the 
loss of heat from the interior. Sources of extra at- 
mospheric heating in c lude p urely hydrodynamic ones 
(Sh owman fc Guillotl 120021 ) and magnetohydrody- 
namic (MHD) ohmic heating due to atm ospheric cur- 
rents in a partially ionized atmosphere 



erna 



et al.l 



l2lfi(jalE iMen^ The 
latter process might not only influence the evolu- 
tionary cooling of planets, but also help govern 
their present-day weather patterns and wind speeds, 
and was suggested simultaneously with, and inde- 
pendently of, t he similar mechanism in the interior 
introduced by IBatvgin fc Stevenson! . Tides might 
also provide another source of atmospheric heat- 
ing, since it is currently unknown whether the dis- 
sipation of tides can deposit a significant amount 
of power above the radiative-convective boundary, 
which deline ates, in our p arlance , "interior" from "at- 
mosphere" dLubow et all 119971; lOgilvie fc Linl 120071; 
IGoodman fc Lacknerll2009f h 

In this paper, we examine the generic character 
of a variety of potential mechanisms that have been 
suggested to explain the inflated hot Jupiters. In ij2j 
we briefly review observed trends in the dependence of 
planetary radii and bulk entropy on mass and incident 
stellar irradiation. In $31 we discuss how atmospheric 
heating affects evolutionary cooling and, hence, radius. 
In Sj4l we discuss how consistently coupling the day and 
night sides of planets in evolutionary cooling calcula- 
tions affects the predicted radii given an extra heating 
luminosity, or the required extra luminosity to match 
a given radius at a given age. In we explore the 
effect of atmospheric ohmic heating, which might be 
a particularly important mechanism in planets with 
large-scale magnetic fields and fast winds. We present 
some details of model atmospheres of a variety of hot 
Jupiters, including free-electron fraction and conduc- 
tivity, and examine how these atmospheric model de- 
tails might depend on gross properties of the planets, 
such as surface gravity and incident irradiation. This 
section concludes by evaluating the stability of planets 
against a potential runaway inflation process. Finally, 
we summarize our conclusions in $6] 

2. DEPENDENCE OF RADIUS AND ENTROPY ON MASS 
AND STELLAR FLUX 

It is instructive to examine the relationship be- 
tween planetary radius and various potentially ex- 
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Fig. 1. — Radius vs. incident flux. For each transiting planet in the radius range from 0.9 to 2.07 Rj found on http:/ /exoplanet.eu 
(circa February, 2013), we plot the planet's radius against the incident stellar flux at the substellar point; points are colored according 
to planet mass. The planets with the largest radii tend to be highly irradiated and low mass. 



planatory variables, such as planet mass and incident 
stellar flux, among the observed transiting exoplanets. 
Here, we simply examine the data for all transiting 
hot Jupiters, making no attempt to correct for selec- 
tion effects. 

Figure [T] presents a scatter plot of planet radius 
versus incident stellar irradiat ion. For each transit- 
ing planet found on http:/ /exoplanet.eu with radius at 
least 0.9 Rj, we color the point by the planet's mass. 
The radii range from our lower cutoff to more than 
2 Rj, and the incident fluxes (Fq values) span more 
than two decades, from fO 8 erg cm~ 2 s _1 to above 
f0 10 erg cm' 2 s -1 . The cloud of points exhibit broad 
scatter, but several trends are evident. Perhaps the 
most obvious of these is that, even though the color 
scale is logarithmically spaced in planet mass, there 
are still far more blue points than red points, indi- 
cating that lower-mass planets are significantly more 
numerous in the sample than higher-mass planets, 
as has also been se en in some slices of Kepler data 
(jHoward et al.ll2~012f ). Furthermore, the planets with 
the largest radii tend to be highly irradiated and rel- 
atively low mass; there are no planets with masses 
greater than ~5 Mj and radii greater than 1.29 Rj. 
The most massive planets do not have radii much 
above 1.0 Rj. Finally, there is a general upward slope 
associated with the cloud of points, suggesting that F 
could be a contributing f actor to explaining planetary 
radii. iDemorv fe Seagen (|2011| ) noted this trend in Ke- 
pler data, as well, finding that incident stellar flux is 



positively correlated with planet radius at irradiation 
levels above ^2 x 10 s erg cm -2 s _1 . 

Figure [5] shows inferred planetary specific entropy 
per nucleon plotted against planet radius. For the 
same sample of planets shown in Fig. [1] we compute 
the bulk interior specific entropy that corresponds to 
the listed radius and mass, assuming solar composi- 
tion and no heavy-element core. Both panels show 
the same cloud of points; the points are colored in 
the left panel by mass and in the right panel by inci- 
dent stellar flux. In the right panel, iso-mass contours 
are overlaid. The transiting planets appear to occupy 
a fairly narrow swath of the radius-specific entropy 
plane, although the true scatter might be slightly dif- 
ferent because of variations in bulk composition and 
core mass. 

3. THE BASIC EFFECT OF ATMOSPHERIC HEATING 
ON RADIUS 

3.1. Modeling Evolutionary Cooling 

At a given mass, the radius of an isolated H/Hc 
sphere with a specified helium fraction depends on ly 
on its interior entropy (jZapolskv fc Salpeterl [TOG gir! 
At the same mass and specific entropy, the radius 
be smaller if there is a heavy element core 



can 



4 Entropy is a useful diagnostic variable because brown 
dwarfs and giant planets are fully convective and, hence, 
at approximately cons tant specific entropy 
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Fig. 2. — Entro py vs. radius, assuming solar composition and no heavy-element core. For each planet in the radius range from 
0.9 to 2.07 Rj on http://exoplanet.eu (in Feburary, 2013), we infer the bulk interior specific entropy that corresponds to the listed 
radius and mass. Top: t he colors of points indicate planet mass. The most massive planets do not have radii much above 1.0 Rj\ 
equivalently, the most inflated planets tend to have lower masses. Bottom: Points are colored according to the incident stellar flux. 
Contours of constant mass in the radius-specific entropy plane are overlaid as green curves. 
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(jGuillot et alj|2006t IBurrows et aLll2007al) or if either 
the helium fraction is large r or there is nonzero metal- 
licity in the bulk i nterior ([Zapolsky & Salpctcr 1969; 
Spie gel et al.ll20lil ). The thermal state evolves with 
time as energy flows from the interior to regions from 
which the optical depth to infinity is low, at which 
point energy eventually radiates to space, leading to 
a loss of entropy and a shrinking radius. This evolu- 
tion depends on the atmosphere, since a higher opacity 
atmosphere will cause a planet to cool more slowly. 

If a planet is not isolated, but rather is irradiated 
by a nearby star, energy propagates inward as well as 
outward. In such a situation, the thermal evolution of 
the interior depends on the net outward radiative flux 
from the deep interior (i^not) through the atmosphere. 
The dayside of a highly irradiated planet tends to have 
a lower net cooling flux for a given mass and entropy 
than an isolated planet has. The effective temperature 
(T c ff) of an object is the temperature corresponding to 

the bolometric net outward flux (T e fj = {Triet/csB} 1 ^ 4 , 
where ctsb is the Stefan-Boltzmann constant). Note 
that the effective temperature, which characterizes the 
difference between the outgoing and incoming fluxes, 
should not be confused with the "equilibrium temper- 
ature" T cq , which characterizes just the incoming flux 
and is roughly the temperature of the photospheric re- 
gion of the atmosphere. In hot Jupiter atmospheres, 
the approximate flux ratio (T eq /T e fj) can be of or- 
der 10 4 . In our evolutionary calculations, we establish 
a mapping between surface gravity (<?), specific en- 
tropy (S), and effective temperature by precalculating 
a large grid of one-dimensional, non-gray, radiative- 
convective atmosphere models with different effective 
temperatures and surface gravities and finding the in- 
terior adiabat associa ted with each one, as described in 
IBurrows et alj (|1997l ). One commonly used approach 
is to take the downward flux in the atmosphere mod- 
els to be the average (over the planet's surface area) of 
the irradiating flux, which is 1/4 of the incident flux 
at the substellar point (this is sometimes described as 
using a "geometric b eaming factor" / = 1/4 - see, 
e.g., the appendix of Spiegel fc Burrowi l2010f h We 
find the function s[T e g,g] and invert it to obtain the 
cooling rate T e s[s,g]. This mapping of (s,g) pairs to 
T e ff is what we refer to as the "atmospheric bound- 
ary condition." (We describe another approach, which 
couples the cooling through different portions of the 
atmosphere — e.g., the day and night sides — in <j4j) 

3.2. The Influence of Atmospheric Heating: Extra 
Absorber 

The temperature-pressure profile of a planet's at- 
mosphere is inextricably linked to the planet's rate of 
cooling. Radiative processes that change the vertical 
thermal structure, therefore, affect the thermal evolu- 
tion. 

A striking feature that has emerged from stud- 
ies of the emergent radiation from hot Jupiters is 
that many of these objects appear to have ther- 
mal inversions in their upper atmospheres, in which 



the temperature increases with heig h t above a rel- 
ative m inimum dHubenv et alJ l2003t IBurrows et al.l 
2007bl; iKnutson et al l 120081: iFortnev et all l2008ab 
Spiegel et al.l l2009t iMadhusudhan & Seageii [2010; 
Knutson et all 120101 : IMadhusudhan! 120121 ). The ther- 
mal inversions presumably result from an enhanced 
opacity in the short wavelength part of the spectrum 
(optical and ultraviolet), above the alkali metal opac- 
ity (mostly sodium and potassium) that some equi- 
librium chemistry models have suggested might con- 
tribute the b ulk of the optical opacit y in hot Jupiter 
atmospheres (iSharp &: Bur rows 200?]). Titanium ox- 
ide (TiO) has been suggested as the source of the 
extra opacity (|Hubenv et al.l 12003). though it is not 
clear whether hot Jupiter atmospheres are vigorously 
enough mixed to bring enough of a heavy, refrac- 
tory species to the upper atmosph ere where it would 
need to be to cause the inversions ([Spiegel ct al. 2009; 
iParmentier et al.ll2013l ). 

If all else about a model is held constant, an en- 
hanc ed optical opacity ac ts as an anti-greenhouse ef- 
fect (H ubenv et alj I2003T ). heating the upper atmo- 
sphere and cooling the deeper atmosphere^ In par- 
ticular, at fixed T G ff, increasing the optical opacity de- 
creases the temperature at depth, in the vicinity of 
the radiative-conyectiv e boundary (see, e.g., Figs. 1— 
5 of lSpiegel et al.1120091 ). and, therefore, decreases the 
entropy of the convective region. However, if we imag- 
ine taking a planet that has no thermal inversion and 
adding an extra optical absorber so as to create a 
thermal inversion, the interior entropy clearly does 
not suddenly decrease in response to altered opac- 
ity. Instead, in order to keep the interior entropy 
constant, the effective temperature (i.e., the interior 
cooling rate) must increase. One might therefore ex- 
pect that models with an extra optical absorber and a 
thermal inversion would cool faster than those without 
these features. 

Figure [3] displays how one-dimensional models in- 
deed exhibit this behavior. The left panel of Fig. [3] 
shows the influence of adding an extra optical ab- 
sorber. At fixed T e ff, the deep atmosphere is cooled 
relative to an atmosphere without the extra absorber 
(represented as k' in the figure). If the effective tem- 
perature is allowed to vary so as to maintain fixed 
interior specific entropy, the model with the extra ab- 
sorber achieves an effective temperature nearly 20% 
greater (215 K vs. 180 K), implying a net cooling rate 
more than twice as great. This increased cooling rate 
is clear in the right panel of Fig. [3l which shows the 
evolution of radius for 1-D model planets with and 
without an extra absorber. The model with the extra 
absorber shrinks more rapidly0 



5 The greenhouse effect, in contrast, results from enhanced 
infrared opacity and heats the deeper atmosphere while cooling 
the upper atmosphere. 

6 However, the effect of an extra optical absorber almost en- 
tirely disappears in the context of a "l+l"-dimensional model, 
as discussed in ij4] below. 
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3.3. The Influence of Atmospheric Heating: Below 
the Photosphere 

Anything that modifies the radiative processes of 
the atmosphere can affect evolutionary cooling. In ad- 
dition to changes in opacity, as discussed in t j3.2l above. 
there are other ways that atmospheres can change the 
mapping of S and g to T e ff. Adding heat to the at- 
mosphere can affect the interior entropy and, there- 
fore, the radius. Consider if some fraction e of the 
intercepted radiation is converted to another form of 
energy (say, mechanical energy) that then dissipates 
as heat below the optical photosphere. At a given 
mass, effective temperature, and incident irradiation, 
a planet where e is greater will have an atmosphere 
that matches to a higher entropy adiabat in the con- 
vective region, and, therefore, that corresponds to a 
higher entropy (and larger radius). 

Figure 0] illustrates how a power source in the 
radiative part of a planet's atmosphere can modify 
the planet's bulk entropy (and, hence, the radius) 
at a given effective temperature. In this figure, a 
range of model planets is displayed, all of which have 
the surface gravity and incident stellar irradiation 
of HD 209458b, and all have an internal flux tem- 
perature of T c ff = 180 K. The atmosphere models 
are calculated with t he atmosphere radiative trans- 
fer code COOLTLUSTY (IHubenvl Il98l IHubenv fc Land 
IT9951: IHubenv et all 120031: Burrows et al.l 120061) . and 
represent solutions to the radiative transfer equation 
in which we enforce radiative and chemical equilib- 
rium. 

In the model planets shown in Fig. |H a small frac- 
tion of the incident irradiation is deposited at various 
depths in the atmosphere, with the power spread in a 
Gaussian distribution over a region of width 0.5 in log 
pressure. The models differ according to the fraction 
of incident irradiation that is deposited at depth in 
the atmosphere (shown on the abscissa) and accord- 
ing to where in the atmosphere the power is deposited 
(shown as different colored curves) . When the power is 
deposited high in the atmosphere (at pressures below 
~1 bar), it has little influence on the interior adiabat 
and, therefore, little influence on the radius. When 
power is deposited deeper, it can have a significant 
influence on the interior adiabat 's entropy, and a cor- 
respondingly large influence on a planet's radius (at 
fixed T c ff = 180 K, with these effects both being larger 
when more power is deposited). The black dashed line 
shows the radius (wl.25 Rj) that this model planet 
has with no extra atmospheric heating. Depositing 
1% of the incident power at 10 bars means that a 
planet with T e s — 180 K has a radius ~5% larger 
(which corresponds to an increment in specific entropy 
of ^0.3 fcs baryon -1 , where ks is Boltzmann's con- 
stant. If the power is deposited at 30 bars, the planet 
has a radius nearly 20% larger, and if the power is 
deposited at 100 bars, the radius is nearly 1.9 Rj. 

Holding T c ff fixed, however, means that the differ- 
ent models in Fig.|4]would correspond to different ages 
in a true evolutionary calculation, with the larger ra- 



dius models corresponding to younger planets. This 
figure demonstrates the general character of the effect 
of atmospheric heating on radius, but a more useful 
comparison would be to hold constant, not T e ff, but 
rather age. It turns out that in the context of 1-D 
models, an extra power source in the deep atmosphere 
can result in a dramatically reduced cooling rate and, 
hence, significantly inflated radii. This effect is muted 
when night-side cooling is taken into account (see SJU). 



4. EVOLUTIONARY COOLING WITH CONSISTENT 
DAY/NIGHT COUPLING 

4.1. The Influence of Non-Isotropic Heating/Cooling 

A gas giant's evolutionary cooling and shrinking 
are mediated by details of the three-dimenisional 
atmospheric structure. Isolated (or nearly isolated) 
objects, such as wide-separation planets or brown 
dwarfs, may, to a reasonable degree of approximation, 
be treated as one dimensional (i.e., as spherically 
symmetric). Hot Jupiters, on the other hand, are ex- 
pected to be tidally loc ked to their stars (although see 
lArras fe~ Socrates 2010{) and, therefore, to experience a 
strong, permanent asymmetry in irradiation between 
the dayside and the nightside. As described in i j3.ll 
treatments in the literature of the thermal evolution 
of hot Jupiters have often used one-dimensional 
models in which the objects are taken to be irradiated 
by the global average irradiation, (l/4)i*o, where Fq 
is the substellar point flux (the "st e llar constant") 
for the planet (iGuillot et al.l 119961: iBurrows et al.l 
2000; Bodcnhcimer ct al. 2001; Burrows et al. 



2003; 

2004 iBaraffe et al 
Lauehlin et al. 
Burrows et al.l 



Mar le v et al 



2008; 




Baraffe et al.ll2003tlBodenheimer et al.ll2003t Gu et al 



Burrows et al.l 12004 iFortnev fc Hubbard! 



2004 T iChabrier et al l 
Baraffe et al.l I2005L 
IFortnev fc Marlev l 
Cha brier fc Ba raffe 



2007; ... .. 

Baraffe et al.l 120081; llbgui fc Burrows! 



Liu et all ._ _ . _ 

2009 : iMiller et all 12001 iLeconte et all !2oToF 
Ibeui et all I2010L 1201 1D ~ If horizontal variations 



in cooling are large, however, these one-dimensional 
models might not accurately capture the character of 
cooling on a highly-irradiated, tidally-locke d planet 
(IGuillot fc Showmanll2002l: IBudai et al.ll2012H 

Hot Jupiters are presumed to be tidally locked in 
synchronous rotation, with one side in permanent day 
and the other in permanent night. The dayside has 
strong static stability imposed by the steady injection 
of energy (and entropy) high in the atmosphere. The 
radiative-convective boundary (RCB) extends deeper 
into the planet as the object cools and shrinks, and 
by an age of 10 9 yrs or more the statically stable re- 
gion on the dayside can extend quite deeply into the 
atmosphere, with a RCB at pressures of hundr eds of 
bars to a kilobar or more (|Burrows et alll200l . The 
nightside, with no irradiation, might have a RCB at a 
much lower pressure (by several orders of magnitude) , 
and, therefore, cool more similarly to an isolated ob- 
ject. We note that s ome dynamical atmospheri c cir- 
culation models (e.g.. iCooper fc Showman! 120051 ) find 
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Fig. 3. — Atmosphere and evolution models of HD 209458b-like planet, with and without extra optical absorber. Top: These 
atmosphere models all have log 10 g = 3.0 where g is measured in cgs units (simi lar to HD 2094 58b). At constant T e ff (dashed black 
curve and blue curve have T e ff = 180 K), an extra optical absorber (re' — see, e.g., Spiegel ct al. 2009) heats the upper atmosphere and 
cools the lower atmosphere. At constant interior specific entropy (dashed black curve and aqua curve have the same interior adiabat 
specific entropy), the extra optical absorber increases T B g by ~20% and, therefore, doubles the net cooling luminosity. Bottom: These 
planets both have a mass of 0.69 Mj, where Mj = 1.899 X 10 30 g is the mass of Jupiter. Comparing the two radius evolution curves 
shows that, as expected, in a 1-D cooling model, the planet with an extra optical absorber and a thermal inversion (blue curve) cools 
faster and its radius shrinks faster than the planet with a "default" (noninverted) atmosphere (black dashed curve). 
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that day-side and night-side temperatures are fairly 
similar by a pressure of ~10 bars, although the poles 
are significantly cooler. Whether the relatively cool, 
less statically stable region is on the nightside or at the 
poles is not terribly important; either deviation from 
an isotropically cooling planet can cause large changes 
to the cooling rate. 

One may, therefore, significantly improve upon 
the models typically used, by appropriately coupling 
simple models of the dayside and nightside together, 
at consistent interior ent r opy and surface g ravity 
(IGuillot fc Showman! l200l IBudai et al.l \20m . In- 
stead of treating the entire atmosphere as though it 
were uniformly irradiated with the planetary aver- 
age irradiation, one may produce a "1+1-dimensional" 
model. First, compute a grid of day-side models, defin- 
ing a function T^ y [s, g). Then, compute a similar grid 

of night-side models, defining a function r c 'g ght [s, g\. 
The total cooling, then, is the sum of the cooling 
through the dayside and of that through the night- 
side: 

-LeS[S,g\ = , (1) 

and T fj [s,<7] is given by the fourth root of the above 
expression. In equation (fTJ, each function must be 
evaluated at the same s and g: the same s because 
convection homogenizes specific entropy, and the same 
g because the planet is (very nearly) spherical. The 
primary effect of strong irradiation on the dayside is to 
create a deep, nearly isothermal layer, through which 
the net flux is nearly zero, thereby strongly decreasing 
evolutionary cooling (i.e., lowering T^ y ). The night- 
side and poles lack the dayside's strong irradiation, 
might remain fully convective out to much lower pres- 
sure levels (bars to tens of bars), and, thus, might more 
efficiently transport heat from the interior to where it 
may radiate away to space. As a result, at equiva- 
lent s and g, the night-side cooling or polar cooling is 
generally much greater than the day-side cooling and, 
therefore, T cS [s,g} 4 ~ (l/2)T" sht [s, g} 4 . This can be 
a large effect that significantly increases the net cool- 
ing relative to what would be found for a hypothetical 
isotropically irradiated planet. 

The process of horizontal redistribution of day-side 
atmospheric heat in hot Jupiter atmospheres has been 
treated in a number of ways in the li terature (see the 
App endix oflSpiegel fc Burrows 2010 fo r a review; see 
alsolCowan fc Ago]||2011l and IMadhusudhan fe Seagerl 
2009). The most physically motivated 1-D parame- 
terization of redistribution is th e P n formalism intro- 
duced by iBurrows et al.l (|2006f ) . in which a fraction 
P n ( "portion-to- night" ) of day-side irradiation is re- 
moved from the day-side atmosphere and deposited 
on the nightside in a specified pressure interval (typ- 
ically ^0.01 — 0.1 bars). P n plausibly ranges between 
(corresponding to no redistribution) and 0.5 (corre- 
sponding to a fully mixed planet where the nightside 
receives half the total intercepted power). In our mod- 
els, we use P n to represent redistribution, removing 



day-side flux and adding heat to the nightside. 

Figure [5] shows the evolution of the radius of 
HD 209458b-like models with various atmospheric 
boundary conditions. Of the seven evolutionary 
trajectories presented, four represent isotropic (1-D) 
models, and three represent day /night- merged (1+1- 
D) models. Two curves are shown that also appear 
in Fig. 13) the "default" 1-D model, shown as a black 
dashed line; and the re' = 0.2 cm 2 g _1 model with a 
thermally inverted layer created by the extra absorber, 
shown as a blue curve. As previously noted, the effect 
of the extra absorber is to increase the cooling rate 
and reduce the planetary radius at a given age. 

However, when we consider the effect of an extra 
absorber in a 1+1-D model with night-side cooling, 
the influence of nonzero re' nearly disappears. The 
day/night-merged models in Fig. [5] are bounded by 
gray dashed-dotted curves that indicate extreme (and 
unrealizable) atmosphere conditions. These cases con- 
sist of isotropic conditions corresponding to day-side 
only (upper curve) or night-side only (lower curve) 
atmospheric conditions. The night-side only object 
has zero irradiation and is significantly heated only by 
the redistributed heat from the dayside, in this case 
corresponding to 30% of the incident irradiation (the 
heating from below corresponding to the net cooling 
flux of the object pales in comparison with the redis- 
tributed heating from the dayside) . Since real planets 
do not experience isotropic day-side or isotropic night- 
side conditions, the top and bottom curves simply il- 
lustrate bounding cases for the possible range of evolu- 
tionary tracks for combinations of irradiated and non- 
irradiated boundary conditions (with P n = 0.3). The 
red and magenta curves show 1+1-D evolution models 
for P n = 0.3 and P n = 0.1, respectively. The green 
curve represents a 1+1-D model for re' = 0.2 cm 2 g _1 
and P n = 0.3. Importantly, the red, magenta, blue, 
and green curves are nearly indistinguishable after a 
few tens of Myr. In other words, when night-side cool- 
ing is incorporated, all 1+1-D models in Fig. [5] cool 
faster than the default model, and subtleties such as, 
in this context, the presence of an optical extra ab- 
sorber become insignificant. 

4.2. The Influence of Atmospheric vs. Deep Interior 

Heating 

Various mechanisms might deposit at depth an 
amount of power that is small in comparison with 
the irradiating power. This extra heating might occur 
deep in the atmosphere, at pressures of bars to hun- 
dreds of bars, or it might occur in the deep convective 
interior. 

Figure |6] shows radius evolution models of an 
HD 209458b-like planet with and without either deep- 
interior or atmospheric heating. Three models demon- 
strate the influence of extra isotropic atmospheric 
heating set to 1% of the incident irradiation (brown 
curves), in which the extra atmospheric power is de- 
posited in a Gaussian distribution of width 0.5 in log 
pressure, centered on 10 bars (thin curve), 30 bars 
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Fig. 4. — Heating at different pressure levels, at constant T e g. Each of the model planets shown here has the surface gravity and 
incident irradiation of HD 209458b. A small fraction of the incident irradiation (fraction shown by the abscissa) is deposited at 
various depths in the radiative portion of the atmosphere (shown as different colored curves) . The power is deposited over a range of 
pressures of logio width 0.5, centered on the pressure identified in the legend. For instance, the green curve corresponds to planets 
where the power is deposited in a range of ~31 bars centered on 100 bars, and the brown curve to planets where the power is deposited 
in a range of ~3.1 bars centered on 10 bars (so the size of the region where power is deposited is approximately constant in altitude). 
The greater the atmospheric heating, and the greater the depth at which the power is deposited, the higher the entropy of the interior 
adiabat and, therefore, the greater the radius. The dashed black line shows the specific entropy and the radius with no extra power. 
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Fig. 5. — Evolution models of an HD 209458b-like planet with and without night-side cooling. Seven radius evolutionary tracks are 
shown of a planet the mass of HD 209458b (0.69 Mj), five of which are represented with thick lines, corresponding to a 1-D model 
(isotropic irradiation - black dashed curve), a 1+1-D model with P n = 0.3 (red curve), a 1+1-D model with P n = 0.1 (magenta 
curve), a 1-D model with k' = 0.2 cm 2 g — 1 (blue curve), and a 1+1-D model with k 1 = 0.2 cm 2 g — 1 and P„ = 0.3 (green curve). (The 
models described as "D/N Merged" are the 1+1-D models with day-side and night-side boundary conditions merged consistently.) 
For illustrative purposes, we include thin gray dashed-dotted lines that indicate evolutionary trajectories for hypothetical planets 
with boundary conditions corresponding to day-side only and night-side only boundary conditions (i.e., all 4tt steradians are bounded 
by either the dayside or the nightside). Real planets clearly do not experience isotropic day-side or isotropic night-side conditions. 
The top and bottom curves represent extreme cases that bound the possible range of evolutionary tracks for mergers of irradiated 
and non-irradiated boundary conditions (with P n = 0.3). Top: The evolution of the seven models is shown measured in Jupiter radii. 
Bottom: The evolution of the same models is shown in units of the radius of the "default" model of the same age. 

(medium curve), and 100 bars (thick curve). Note 
that when extra heating is at a depth of tens of bars or 
more, it acts to retard evolutionary cooling, in contrast 
to the effect of k' — which is essentially extra heat- 
ing at millibars or higher — which acts to accelerate 
cooling and shrinking. Differences in radius between 
different models generally (but not always) shrink with 
time, but after 1 Gyr, the model in which 1% of in- 
cident irradiation is deposited at 10 bars is ~0.1 Rj 
larger than the default model; the model with extra 
power at 30 bars is more than 0.2 Rj larger than the 
default model after 1 Gyr; and the model with extra 
power at 100 bars is more than 0.5 Rj larger than the 
default model after 1 Gyr, and has a radius larger than 
1.7 Rj at this age. This is consis tent with the find- 
ings of iGuillot fc Showman! ([2002). who argued that 
depositing 1% of the incident irradiation deep in the 
atmosphere could have an important radius-inflating 
effect. 

Models in which the dayside's extra atmospheric 
heating is merged with night-side cooling exhibit a sig- 
nificantly reduced influence of extra day-side heating. 
This coupled day /night cooling scenario with day-side 
heating is represented (for the case of heating centered 
at 100 bars) by the purple curve. If the nightside is 
unheated, the extra heating on the dayside, which re- 



sults in such dramatically increased radii in the 1-D 
cases, has essentially no influence on the radius in the 
1+1-D case (note the very small difference between 
the 1+1-D models with no day-side heating and with 
day-side heating, but no night-side heating — the red 
and purple curves, respectively). 

Figure[6]also presents two examples of radius evolu- 
tion with deep-interior heating — one for a 1-D model 
(green) and one for a 1+1-D model (aqua). The in- 
terior heating in both models is set to 1% of the in- 
cident irradiating power and is assumed to occur in 
the center of the planet. Note that an extra power 
source that is deep in the envelope, but not in the cen- 
ter, might have a complicated influence on the ther- 
mal structure of the object. Instead of resulting in 
uniformly higher specific entropy in the deep interior, 
such an extra power s ource might inhibit convection 
(jWu fe Lithwickl 12013). We avoid this complexity by 
assuming that the extra heat is deposited in the cen- 
ter, but we note that fully exploring the structural 
influence of deep envelope heating is a subject ripe 
for future work. At early times, the brown curves (1- 
D models with atmospheric heating) have larger radii 
than the models with an equivalent amount of cen- 
tral heating. However, the models with deep-interior 
heating reach an inflated asymptotic radius because 
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~0.1 Rj, such that, in the set of models displayed in 
Fig. achieving the same asymptotic radius requires 
twice as much central heating power in the context of 
a 1+1-D model with P n = 0.3 as with a 1-D isotropic 
model. 

4.3. The Influence of Redistributive Winds (P n ) 



Fig. 6. — Evolution models of an HD 209458b-like planet, with and without deep-interior or atmospheric heating. Eight models 
are presented: the default (isotropic heating) model (black dashed); three 1-D models with extra atmospheric heating set to 1% of 
the incident irradiation, deposited at 10 bars (thin brown), 30 bars (medium brown), and 100 bars (thick brown); a 1+1-D model 
with P n = 0.3 and with extra heating deposited at 100 bars and an unheated nightside (thick purple); a 1+1-D model with no extra 
heating and P n = 0.3 (red); a 1-D model with deep-interior heating that is 1% of incident (dark green); and a 1+1-D model with 
central heating that is 1% of incident and with P n = 0.3 (aqua). In order to limit the number of curves for ease of viewing, we 
present only a single purple curve (i.e., we do not display the evolution for when the heating is deposited at different pressure levels). 
Depositing the heat at 10 or 30 bars makes no discernible change to the purple curve. With 1-D models, if extra heating is deposited 
deep (e.g., 100 bars), the radius can be inflated more by heating in the atmosphere than by heating in the deep interior until late 
times (several Gyr); however, when night-side cooling is incorporated, extra heating in the deep interior inflates the radius more at 
all times later than a few Myr. 

in such a model the effective temperatures has a strict 
(positive-definite) floor corresponding to the flux im- 
plied by the power source. In contrast, the asymp- 
totic radius of a model with extra atmospheric heat- 
ing is the zero-temperature, fully degenerate radius — 
approximately 1 Rj for objects of roughly Jupiter's 
mass — because the extra heating simply reduces the 
net cooling rate to a lower, but still nonzero, value. 
The asymptotic radius of the 1-D model with central 
heating is roughly ^0.08 Rj larger than that of the 
1+1-D model with central heating and P n = 0.3. 

Figure [7] depicts the radius evolution of models of 
an HD 209458b-like planet with a range of levels of 
central heating power. This figure particularly high- 
lights the contrast between 1-D models (green) and 
1+1-D models with P n — 0.3 (aqua). For both the 
isotropic model and the one with night-side cooling, in- 
creasing the interior heating power (measured in units 
of the incident power) from 0.5% (thin curve) to 1% 
(medium curve) results in nearly a 0.1 Rj increase 
in asymptotic radius, and increasing the power from 
1% to 2% (thick curve) results in another ~0.1 Rj 
increase. For this set of models, including night- 
side cooling with P n = 0.3 results in a decrement of 



As described in §4.11 in the absence of winds that 
redistribute energy from the dayside of a planet to its 
nightside, the nightside would be much cooler than 
the dayside. Day-night temperature differences, how- 
ever, lead to pressure gradien ts that drive redist ribu- 
tive winds (see, e.g., Fig. 3 of lBurrows et al.l feOlO'). 

Figure [8] demonstrates both the effect of incorpo- 
rating consistent night-side cooling in an atmosphere 
model and the influence of the P n redistribution pa- 
rameter. For variety, in this figure we show the evo- 
lution of a model with the mass and irradiation of 
WASP-12b (IHebb et alJl2009h instead of HD 209458b. 
The model planet is evolved with three different atmo- 
spheric boundary conditions and with the extra inte- 
rior p ower that is necessary, according to llbgui et al.l 
(2010), to maintain its current radius at Rj 
(2.4 x 1O~ 6 £ , or 0.2% of the incident power). As 
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Fig. 7. — Radius evolution models of an HD 209458b-like planet with deep-interior heating, with and without night-side cooling. 
Seven models arc presented: the default model, and six models with central heating — 0.5%, 1%, and 2% of incident (thin, medium, 
and thick curves, respectively). Three of the models with central heating have isotropic boundary conditions (dark green), and three 
have 1+1-D boundary conditions with P n = 0.3. In this set of models, including night-side cooling with P n = 0.3 means that the 
model requires a factor of ~2 more interior heating to maintain the same radius. 



above, the source of this extra interior power is not 
specified in this model, but it could correspond to 
tidal heating or to ohmic heating in the interior. The 
three boundary conditions correspond to a default 
(isotropic) 1-D model (/ = 1/4), and two 1+1-D 
models (with redistribution given by P n = 0.1 and 
P n = 0.3) that couple day-side cooling to night-side 
cooling in accordance with Eq. (JTJ). Redistribution of 
heat from the day to the night has the dual effects of 
cooling the dayside and heating the nightside. The for- 
mer effect slightly increases the dayside's cooling rate, 
but the latter effect significantly slows night-side cool- 
ing. The influences of these two effects on radius evo- 
lution have opposite sign, but the (radius-boosting) ef- 
fect of a warmer nightside is quantitatively much more 
significant than the (radius-shrinking) effect of a cooler 
dayside. As a result, the 1+1-D models with consis- 
tent day /night cooling cool more rapidly than does 
the 1-D isotropic model, and the model with P n — 0.3 
asymptotes to a radius more than 0.1 Rj smaller than 
the model with isotropic conditions Q The isotropic 
model, unsurprisingly, asympt otes at the same radius 
as the analogous model from llbgui et all ()2010f ) with 
identical interior luminosity. The model with P n — 0.3 
has a warmer nightside than the one with P n = 0.1, 
which leads to a lower net flux (i.e., a lower effective 
temperature). Therefore, since the model planet's net 

7 This same effect was seen previously in Figs. \E\— for 
HD 209458b. 



cooling is dominated by the nightside's effective tem- 
perature, the P n = 0.3 model asymptotes at a larger 
radius (by 0.03 Rj) than the P n = 0.1 model. 

5. THE EFFECT OF OHMIC HEATING 

5.1. Heating in the Atmosphere 

The temperatures and densities that are present 
in the day-side atmospheres of highly irradiated 
planets dFortnev et al.l l2008at ISpiegel et all 120091: 
Mad husudhan fc Seagerl 120101 : ISpiegel fc Burrowsl 



2010|j imply a small, but nonnegligible, free electron 
fraction from the partial ionization of species with 
low ionization potential (mostly atomic potassium 
and sodium). During the last decade, a wide variety 
of circulation models have predicted that, near the 
equator and at depths comparable to the optical and 
infrared photospheres, hot Jupiters should have east- 
ward (prograde) wind velocities (u) of ~1 km s -1 or 
more dShowman fc Guillot|l2002T iCooper fc Showman! 
1 20051: jLangton fc Laughlinl 120081: Tbobbs-Dixon fc Linl 
2008 : iMenou fc Rauscherll2009t iShowman et al.l l2009: 
iHeng et al.1 120111 ). If planets have large-scale dipolar 
magnetic fields (B), zonal winds will carry the 
free electrons across field lines, thereby inducing 
an electric field in respon se to the no n zero u x B, 
as pointed out by b oth iPerna et al.l ()2010aD and 
iBatvein fc Stevenson! (|2010f ). The electric field will 
drive curren t loops that could close either in the 
atmosphere (jPerna eTaTl I2010al) or in the interior 
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Fig. 8. — Radius of WASP-12b-like planet with interior heating. The evolution of radius is shown for a model planet the mass of 
WASP-12b (1.39 Mj) for three different boundary conditions: a default (/ = 1/4) model, and two models where the day and night 
sides are consistently coupled, with redistribution parameter P n set to 0.1 and to 0.3. In all three cases, there is an extra interior 
power source equal to 2.4xl0 — 6 L© (the value needed to maintain WASP-12b's radius, according to Ibgui ct al. 2010). The planet 
cools more through the nightside than through the dayside (i.e., the effective temperature of the nightside is greater at equivalent 
entropy and surface gravity), which causes the models that incorporate ni ght-side cooling to reach a lower asymptotic radius than 
is reached in the default / = 1/4 case (which is essentially the model of Ibgui ct al. 2010). Greater redistribution increases the 
night-side atmospheric temperatures, which slows the night-side cooling. Since the cooling is dominated by the night-side cooling, 
this also reduces the total cooling rate and leads to a slightly larger asymptotic radius at equivalent central heating power. 



iBatygin fc Stevenson! (|2010D . Here, we consider those 
that close in the atmosphere. 

Figure [9] shows five day-side atmospheric profiles 
(at fixed T^ff) for a planet with the incident stellar ir- 
radiation of HD 209458b. One profile is for a planet 
with no extra heat ing; the other four have heating 
profiles taken from iPerna et al.l (|2010b[ ). for 3- Gauss 
(G) magnetic fields (with and without associated drag 
on the flow) and for 10-G magnetic fields (with and 
without associated drag) . Regions of the atmosphere 
that are co n yective are indicated in green. In the 
IPerna et al.l (|2010b| ) models, stronger magnetic (and 
zero drag) fields lead to greater atmospheric heating. 
As seen in Fig. [U greater heating (at fixed T c g) leads 
to greater interior entropy and larger radii. Figure[9]il- 
lustrates how atmospheric heating at constant T e g can 
modify the temperature-pressure profile to cause the 
atmosphere to match a higher entropy adiabat — the 
same process that leads to the radius inflation effect 
of atmospheric heating seen in Fig. [6] 

5.2. Constraints on Ohmic Dissipation 

The ultimate power source of ohmic dissipation is 
stellar irradiation. Heating from the star drives winds, 
which, with a planetary-scale dipolar magnetic field 
and with free electrons in the atmosphere, lead to cur- 
rents in the atmosphere that ohmically dissipate, drag- 




1500 2000 2500 3000 3500 4000 4500 
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Fig. 9. — Temperature-pressure profiles in the presence 
of ohmic heating. Five atmospheric profiles are shown (at 
fixed T c ff) for a planet with the incident irradiation of 
HD 2094 58b. We ad o pt the atmospheric heating profiles pre- 
sented in[Pcrna ct al. (20103) for 3-G magnetic fields (with and 
without associated drag on the flow) and for 10-G magnetic 
fields (with and without associated drag). Stron ger magnetic 
fields l ead to greater atmospheric heating in the IPerna et al.l 
(2010b) models, and greater heating causes an atmosphere with 
a given effective temperature to have a base adiabat that is at 
higher entropy. Green regions are convective (and have uniform 
entropy) and other regions are radiative (and have positive ra- 
dial entropy gradient). 

ging the winds and releasing heat in the atmosphere 
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or in the convective interior. The power released by 
this process is bounded by the power incident on the 
planet. 

The ohmic heating power per unit volume is 
dPohm/dV = J 2 /cr, where J is the current density and 
a is the electrical conductivity. For a dipolar magnetic 
field stren gth B, with zo nal wind speed it, this can be 
expressed (Menou 2012j) as 



dp 



ohm 



dV 



(2) 



where c is the speed of light. The electrical conductiv- 
ity can be written 



n e e 2 /m e 
n n (Av) 



(3) 



where, n e and n n are the number densities of free 
electrons and neutrals, respectively, e is the electron 
charge, m e is the elec tron mass, A is th e electron 
cross section, and, per iDraine et al.l (|l983), {Av) ~ 
10" 15 (128fc B T/97rm e ) 1 / 2 . The ohmic heating density, 
therefore, is 



dP ( 



ohm 



-3 „-l 



dV 



• 10 erg cm s 



B 



lkm/s/ V 10G 
Ye \ ( T x ~ 1/2 



10- 4 J \10 3 K 



(4) 



where the free-electron mixing ratio Y e = n e /n n . The 
scaling for Y e (10 -4 ) might seem strangely high, since 
at common conditions in hot Jupiter atmospheres 
most free electrons are contributed by singly ionized 
potassium and sodium, whose abundances (mixing ra- 
tios of ~10 -7 and ^2 x 10~ 6 , respectively, if the com- 
position is solar) are both much lower than 10 -4 . For 
the hottest hot Jupiters, however (e.g., WASP-12b and 
HAT-P-7b), ionization of hydrogen can contribute up 
to ~10~ 4 in free-electron mixing ratio (see Fig. [TU|) . 

If the ohmic heating occurs over a solid angle it 
of the planet's dayside, and extends vertically over N 
scale heights H, then the integrated ohmic heating 
rate is 



(5 x 10 30 ^ 
V s 



Rj 



B 



H 



350 km 



lkm/s/ V 10G / V 10 ~ 



T 



10 3 K 



(5) 



-1/2 



The power incident on a planet a distance a from a 
star of luminosity i* is 



P in ~ 7 x 10 30 erg s" 1 



(o.Ol Au) 



( 6 ) 

Since the ohmic heating clearly cannot exceed the in- 
cident stellar heating, a comparison of equations ([5]) 



and (0 shows that there must be feedback between B 
and u, whereby ohmic dissipation limits wind speeds. 
Since the heating power comes from the d issipation of 
wind v el ocities, this point is trivially true (|Perna et al.l 
l2010al rbl: lBatvgin et allfeoill ). Nevertheless, it is inter- 
esting to note that simple energy balance requires that, 
if the dipolar planetary magnetic field is of order 10 G, 
the hottest hot Jupiters cannot (over a large fraction 
of their day-side atmospheres) have zonal winds at the 
^kilometer-per-second level that is typically suggested 
by numerical simulations that do not include the mag- 
netic drag effect. 

[Menou (2012) suggests scaling relations for winds 
in ohmically-heated hot- Jupiter atmospheres. With- 
out considering ohmic dissipation, the horizontal mo- 
mentum equation leads to an o rder-of-magnitude esti- 
mate of ((Showman et al. 2010h 



fcAT, 



horiz 



firrip 



AlnP. 



(7) 



where AXhoriz is the day-night temperature contrast, 
fi is the mean molecular weight of the atmosphere 
in units of the proton mass m p , and AlnP is the 
number of scale heights above a base level of the at- 
mosphere at which the day-side and night-side atmo- 
spheres have the same thermodynamic state. Note 
that equation ((7j) may be expressed as 



( C sday 



5 night 



AlnP. 



"ff night l 



[2-kH]) AlnP, (8) 



where the c s terms are the sound speeds on the 
day and night sides, and the c g {2irH] terms are 
the gravity wave speeds for wavelength 2irH on the 
day and night sides. For reasonable values of the 
day-night temperature difference, equations (0 and 
((8|) yield zonal wind speeds of order the day-side 
sound speed, roughly a few kilometers per second 

(~2{T day /1000 K} 1/2 km s" 1 ), which indicat es that 
the m agnetic drag term in equation (5) of iMenoul 
(|2012l ) plays an important role in limiting wind speeds 
if the magnetic field is of order several Gauss or 
more (|Perna et al.ll2010alfbh . Note that most re d-giant 
hot Jupiters (jSpiegel fc Madhusudhanl I2012D prob- 
ably do not receive enough irradiation to have a 
high enough free-electron fraction to experience non- 
negligible ohmic heating. 

5.3. Ionization and Heating in the Population of 
Known Hot Jupiters 

It is now possible to construct radiative-convective 
atmosphere models, guided by observations, of more 
than 20 transiting exoplanets. We have pub- 
lished models of over a dozen, including, in or- 
der of increasing inc ident stellar flux , HP 189733b 
Burrows et all 120081) HP 209458b (IKnutson et al. 
20081 : iBurrows et all 120011. TrES-3 (iFressin et al 



2010 ), HAT-P-7b (ISniegel fc Burrows! I2010H. WASP 



18b (Machalek et aLll2010 submittedl )7and WASP-12b 
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(jCowan et al.H2012l) . The radii of the members of this 
sextet, respectively, are 1.14 Rj, 1.35 Rj, 1.34 Rj, 
1.36 Rj, 1.11 Rj, and 1.83 Rj. 

Figure [10] provides illustrative examples of the 
types of atmospheric structures that may be found 
across a range of hot Jupiter conditions. We present 
vertical profiles of (clockwise from upper left) atmo- 
spheric temperature, electron mixing ratio, conductiv- 
ity, and ohmic heating for the six planets mentioned 
above, for which the incident flux varies from 0.48 to 
9.10xl0 9 erg cm~ 2 s^ 1 , and logio[ff] (in cgs) ranges 
from 2.96 to 4.42. The profiles assume day-side aver- 
age conditions. 

The top-left panel of Fig. [TU] shows that three of 
these models - HD 209458b, HAT-P-7b, and WASP- 
12b - have thermal inversions in their upper atmo- 
spheres, alth ough that of WASP-12b is somewhat 
dispu ted (see IMadhusudhan et al.ll201lL ICowan et al.l 
I2012L and lCrossfield et al.ll2013l for further discussion). 
All one-dimensional models of sufficiently highly irra- 
diated planets have deep, nearly isothermal, layers, 
extending more than two decades in pressure down 
from the photosphere at ~1 bar. In general, the 
temperature of the isothermal layer scales with the 
strength of the irradiation roughly in proportion to 

Fq 1 ^ 4 , where Fp is the inc ident flux at the substellar 
point (|Hubenv et al.l2003h . However, the presence of a 
thermal inversion tends to heat the upper atmosphere 
and cool the lower atmosphere relative to the profile 
that would obtain without an inversion, as is clear 
from examining the relative temperature profiles of, 
e.g., WASP-18b and WASP-12b. The electron mixing 
ratio is calculated via Saha equations, assuming so- 
lar composition of the elements. The conduc t ivity is 
calculated with equation (|3J), found in IMenoul (|2012f ). 

In order to produce the heating profile plot in the 
lower right panel of Fig. [TU1 we use equation Q to 
calculate the ohmic heating density, and vertically in- 
tegrate this density times ttR 2 , from the outside of 
the planet inward. The abscissa, then, is the ratio 
of the cumulative ohmic power (from the outside of 
the planet down to a given ordinate) to the incident 
stellar power. Note that the ratio of the full integral 
of the ohmic power to the incident irradiation is the 
efficiency factor e referred to in 35.41 below. The tem- 
perature and electron mixing ratio profiles are calcu- 
lated from the model atmospheres, and the zonal wind 
speed u and magnetic field strength B are fixed at 
values motivated by previous studies in the literature. 
We present two heating profiles for each planet. The 
solid curve assumes zonal wind speeds of 1 km s _1 
throughout the whole atmosphere, where this speed 
is motivated by the results of a variety of circulation 
models in the literature and by the simple scalings in 
equations and {5} . The dashed curve assumes that 
the zonal wind speed is prograde (relative to the co- 
rotating frame) by 1 km s _1 at pressures lower than 
10 -4 bars, stationary (0 km s _1 ) at pressures greater 
than 300 bars, and varies linearly with altitude in be- 
tween. The magnetic fields are assumed to be 3 G 



for all planets. Note that increasing B from 3 G to 
10 G increases the heating by roughly an order of mag- 
nitude, which results in the cumulative ohmic power 
being greater than the incident stellar power for the 
models of HAT-P-7b and WASP-12b if the winds are 
~1 km s _1 over a significant portion of the day-side 
atmosphere. In truth, the wind speeds would depend 
on the magnetic field strength, which itself depends 
on hydrodynamic motions in the deep interior and 
in the atmosphere that, presumably, lead to a mag- 
netic dynamo. Numerous works have attempted to 
simulate or calculate aspects of this feedback, with 
varying levels of sophistication, a lthough none of the se 
works capture s the full process dPerna et al.l l2010allbl; 
Bat ygin et all 120111: lMenoull2012t iRauscher fc Menoul 
l2013t IWu fc Lithwickl 120131) . The profiles presented 
here should be taken not as quantitative predictions, 
but rather as qualitative guides to the vertical varia- 
tion in heating that might obtain in hot Jupiter atmo- 
spheres at a range of levels of irradiation. These pro- 
files may be scaled with u 2 B 2 for different assumptions 
of wind speeds and magnetic field strengths. Clearly, 
the largest contribution to the degree of ohmic heat- 
ing is the local conductivity, which varies by ~8 orders 
of magnitude across different altitudes and different 
planets. 

The ionization depends very sensitively on temper- 
ature. On the nightsides of these planets (calculated, 
but not shown) , the temperatures are cool enough that 
electron mixing ratios and, therefore, electrical con- 
ductivities, are much lower than on the daysides. The 
night-side conductivities are all < 10 4 s _1 , render- 
ing any implied heating rate too low (<10 -4 of the 
incident stellar power) to have a significant effect on 
planet radius. 

Our one-dimensional atmosphere models of the day 
and night sides of highly irradiated gas-giant plan- 
ets suggest that the nightsides might have very low 
free-electron fractions and experience negligible ohmic 
heating. The daysides of the hottest hot Jupiters, 
however, should be highly enough ionized that that 
they might experience significant ohmic heating if 
wind speeds are as high as dynamical models have 
suggested and if the magnetic field strengths are 
as great as Jupiter's ^ 4-Gauss field (jRusselll 1993; 
iChristensen et al.ll2009f ). 

5.4. Planetary Radii Are Generally Stable Against 
Runaway Expansion Due to Ohmic Heating 

In order to investigate the stability of a planet's 
radius, we consider how an equilibrium situation re- 
sponds to a perturbation. Define R® as the equilib- 
rium radius. That is, at R p = R®, heating balances 
cooling and dR p /dt = 0. Let s be the specific entropy, 
i.e., the total entropy S per unit mass M p : S = J sdm 
and dS = sdm. In a planet's convective interior, s is 
nearly constant. 

Note that a change in heat corresponds to a change 
in entropy vi&dQ — TdS = Tsdm, where the operator 
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Fig. 10. — Profiles of properties of day-side atmospheres of transiting exoplanets. Each panel shows pressure as the ordinate 
(decreasing upward on a logarithmic scale; i.e., altitude increases upward approximately linearly) and some other atmospheric 
property as the abscissa. The profiles are color-sorted by incident flux, and the legend in each panel indicates incident flux ("F:") 
in units of 10 9 erg cm~ 2 s — 1 and the base-10 logarithm of surface gravity ("logg:") in dynes per gram. Top left: Temperature 
profiles are shown, including some corresponding to thermally-inverted atmospheres (e.g., HD 209458b, WASP-18b, and WASP-12b 



erg cm 
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although see [Ivladhusudhan et al. 2011). More highly irradiated planets tend to have hotter isothermal layers between ~1 and 



^300 bars. Top right: Electron mixing ratios are derived from pressure and temperature according to Saha equations. Most of 
the free electrons in the upper atmosphere come from singly ionized sodium and potassium. All atmospheres shown are at most 
weakly ionized (the maximum ionization shown is ~10~ 4 , for a model of WASP-12b), but, importantly, the ionization is nonzero. 
More highly irradiated planets have greater partial ionization. (High in the atmosphere — at ~microbar pressures — there can be a 
nonnegligible photoionization component to the free electron fraction.) Bottom left: Conductivity profiles are derived from the free 
electron number density. Bottom right: Heating profiles are derived from the conductivity, with assumptions about the magnetic 
field strength (3 Gauss) and wind speeds (solid: 1.0 km/s; dashed: 1.0 km/s at pressures lower than 10~ 4 bars, 0.0 km/s at pressures 
greater than 300 bars, and varying linearly with altitude in between). 



d indicates an incompete differential. So, 

Tsdm I dt = dQ = L nc t^ ln dt — —Ldt , (9) 

Mass / 

where L nc t_i n = —L is the net inward (heating) power, 
which is the negative of L, the net outward (cooling) 
luminosity. The net luminosity L = C — H , where C 
is the outward luminosity and H is the inward lumi- 
nosity. 

These are the ingredients needed to evaluate 
dRp/dt — i.e., {dR p / ds){ds / dt) — for a given per- 



Hclmholtz time: 



tkh 



J Tdm 



E 
L 



(11) 



turbation Ai? p from R®. If we define 6r p 
we find that, unsurprisingly, 



dSr. 



dt 



AR P / R p , 



(10) 



where the relaxation timescale tkh is the Kelvin- 



Here, C is a simple function of planet radius and T e ff, 
the effective temperature. H is the surface integral of 
e F, where F is the irr adiating flux, and the efficiency 
s ([Batygin et al.l 1201 If ) with which irradiation is con- 
verted to heating may be thought of as a function of 
R p , a possible magnetic field B, planetary mass M p , 
the irradiation itself (F), and perhaps other variables. 

Finally, E is the product of the number of particles 
(of order ^Mp/m p ) with the mass-weighted average 
of kT. For simplicity, we will here consider e to be a 
function of at most R p (and perhaps constant), while 
keeping in mind that this is a simplification. Again, 
dRp/dt equals zero at R = R because eF-4aT^ s = 
for R p = R^. The form of tkh makes clear that it is 
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essentially the Kelvin-Helmholtz timescale, which is 
the natural timescale on which an object responds to 
thermal perturbations. Importantly, a perturbation 
relaxes to equilibrium since the coefficient of Sr p in 

the relation 5r p cx 5r p is negative. 

Figure [11] shows illustrative curves for heating 
power versus cooling power for models of an irradi- 
ated hot Jupiter planet (similar to HD 209458b). In- 
tersections between the heating curves and the cooling 
curves represent equilibrium radii ("i?p") where the 
cooling balances the heating. The upshot of the sta- 
bility argument just presented is evident in this figure, 
because at all equilibria the slope of the cooling curve 
is greater than the slope of the heating curve, indi- 
cating that a small positive perturbation to the radius 
will lead to cooling exceeding heating and, therefore, 
to the radius shrinking back to the equilibrium. This 
indicates that, in these models, there is no purely ther- 
mal runaway instability in radius (for s ^ 10%). This 
result is in contrast to the suggestion in lBatvgin et al.1 
(|2011[ ) that ohmic heating can lead to a thermal in- 
stability that culminates in complete evaporation of a 
planet. When a planet's radius becomes large enough, 
however, a runaway Roche-lobe overflow process can 
still occur that might initially be triggered by ohmic 
heating. 



6. CONCLUSIONS 

We have presented an analysis of various sub- 
tle, but important, physical effects that can influ- 
ence the radii of hot Jupiters. To explain the anoma- 
lously large radii of an interesting subset of this class, 
some have invoked the deposition of thermal power ei- 
ther in planets' deep, convective interiors or at var- 
ious levels in their radiative atmospheres. To ex- 
plore the consequences of such heating on the radii 
and radius evolution of hot Jupiters, we generated 
self-consistent atmosphere-planet evolutionary mod- 
els and found that the evolution of highly irradiated 
gas-giant planets depends sensitively on where the en- 
ergy is deposited. If "extra" power goes into heating 
the upper atmosphere, via an extra thermal-inversion- 
causing optical absorber, this (slightly) increases the 
evolutionary cooling rate and results in slightly smaller 
radii. If the additional power is deposited more deeply 
in the radiative atmosphere (e.g., via ohmic heating), 
the interior cooling rate decreases, with the result that 
the planet's radius is larger at a given age. The radius- 
boosting effect of a given amount of power increases 
with the depth at which that power is delivered in the 



atmosphere. However, the asymptotic radius in the 
limit of infinite age is not changed. In other words, 
atmospheric heating might be able to produce a large 
inflating effect, but one that does not last forever. In 
contrast, the same power deposited deep in the convec- 
tive interior of the planet always has a larger asymp- 
totic radius than if the same power is delivered any- 
where in the radiative zone. 

Importantly, cooling through the (cooler) nightside 
of a tidally locked hot Jupiter, or through its (cooler) 
poles, increases the rate of interior cooling and radius 
shrinkage relative to isotropic models. As a result, 
the radius-inflating influence of atmospheric heating 
is reduced in a 1+1-D model that consistently treats 
day/night redistribution of heat and cooling through 
the nightside. This makes it difficult to significantly 
increase the radius of a hot Jupiter via atmospheric 
heating alone. 

In addition, we find a limit on the prodcut of zonal 
wind speed and magnetic field strength in hot Jupiter 
atmospheres that have significant ionization. Specifi- 
cally, the most highly irradiated planets cannot have 
uB > 10 km s" 1 Gauss over a large fraction of their 
daysides, where u is the zonal wind speed and B is the 
dipolar magnetic field strength in the atmosphere. 

The vexing problem of the inflated radii of some 
hot Jupiters has persisted for more than a decade. 
It has long been appreciated that intense irradiation 
can slow the evolutionary cooling of these objects, 
whatever other processes are invoked to explain the 
anomaly. Here, however, we have shown that coupling 
the day-side cooling with that through the nightside or 
the polar regions largely undoes the reduced interior 
cooling effect of the irradiating flux. An extra power 
source may, then, be needed in their deep interiors, not 
merely in their atmospheres, to explain most inflated, 
multi-billion- year-old planets. Moreover, the degree of 
extra power that must be deposited in the deep inte- 
rior would seem to be greater than has been inferred in 
those many studies using isotropically cooling models. 
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